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Abstract. We present a systematic classification of the elementary excitations 
of Bosc-Einstein condensates in cylindrical traps in terms of their shapes. The 
classification generalizes the concept of families of excitations first identified by 
Hutchinson and Zaremba (1998) Phys. Rev. A 57 1280 by introducing a second 
classification number that allows all possible modes to be assigned to a family. We 
relate the energy ordering of the modes to their family classification, and provide a 
simple model which explains the relationship. 



1. Introduction 

Collective modes provide important signatures of the physics of Bose-Einstein 
condensates, and they have been intensively studied over the past few years. A large 
number of theoretical calculations, both numerical (e.g. 0]-@) and analytical (e.g. JF|- 
have been made based on the standard Bogoliubov treatment, and experimental 



observations [|L9[|-P4|I have confirmed that at very cold temperatures the Bogoliubov- 
de Gennes (BdG) equations provide a very accurate description of these modes. This 
method amounts essentially to specifying eigenmodes for a perturbative, linearized set 
of equations around a potential formed by the trap and the ground state condensate 
wave function. In this paper we use the standard BdG equations to make a systematic 
study of the shape of elementary modes and the relationship to their energy ordering. 
We note that very recent studies have incorporated finite temperature effects in order 
to improve accuracy at temperatures closer to the critical temperature (e.g. see p5 



28] and references therein), however, for our purposes, the BdG equations provide an 
appropriate and tractable formalism. 

For an isotropic trap, the spherical symmetry allows a separation into a radial 
wave function and spherical harmonics, thus much of the behaviour can be predicted 
from familiar cases in linear quantum mechanics. For axi-symmetric traps, which is 
the case for most experimental situations, such a separation is not possible, and a 
variety of techniques have been used to obtain the mode functions. In an elegant paper, 
Fliesser et al Jl6| obtained approximate analytic solutions in the hydrodynamic limit by 



identifying a non-trivial operator which commutes with both the eigenvalue operator of 
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the BdG equations, and the angular momentum component along the symmetry axis. 
Hutchinson and Zaremba |J obtained numerical solutions for the low-lying modes of 
axially symmetric traps, and studied the behaviour of these modes as the asymmetry of 
the trap was changed from prolate to oblate. The latter authors noticed that the lowest 
lying modes could be grouped into four distinct families, based on the behaviour of the 
mode frequency as a function of the trap asymmetry. They also showed by example 
how some of their modes could be identified with the modes found by Fliesser et al. 

In this paper, we extend the results of Hutchinson and Zaremba, and Fliesser et 
al by obtaining a classification of the modes in terms of their geometrical shapes. This 
enables us to give a generalized definition of the mode families (which are in principle 
of unlimited number). Our classification into families can be expressed in terms of 
well defined quantum numbers, which can be directly related to the quantum numbers 
used by Fliesser et al. In addition, we have examined the energy ordering of the 
modes for the case of oblate and prolate traps, and give a simple model that relates 
the characteristic quasiparticle shapes to their energy ordering. This also explains the 
results of Hutchinson and Zaremba, that all modes in a family have the same frequency 
dependence on trap anisotropy. 



2. Formalism 



The ground state of a dilute Bose-Einstein condensate at low temperatures is well 
described by the Gross-Pitaevskii equation (GPE) 

~H + NU \i,{r)\ 2 ]i,{r) = fi^r), (1) 

where 

H = ~V 2 + V t[ap (r) (2) 
2m 

is the single-particle Hamiltonian of the trap and N the number of atoms. The effective 
interaction strength Uq is related to the s-wave scattering length a and the atomic mass 
m by Uq = 4nh 2 a/m (e.g. see [|L4|). The wave function ip(r) is normalized to unity 
according to J \ip(r)\ 2 dr = 1 and the eigenvalue \i is the chemical potential of the 
condensate. The quasiparticle excitations on the ground state ipg{ r ) have amplitudes 
Ui and Vi determined by the Bogoliubov-de Gennes (BdG) equations 

Cui(r) + U ^l(r)vi(r) = hu^r) 
Cvi(r) + U ip* 2 (r)ui(r) = -hu^r), 

where C = H + 2NU Q \i/j g (r)\ 2 — fx We note that the solutions itj and of equations 
(0) are not necessarily orthogonal to the ground state, although this is required in a fully 
quantum mechanical treatment. However, since the orthogonal excitations can readily 
be obtained from the solutions of equation (§) by projection we will present only 
the direct solutions of equations (H) here. 
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In this paper we will consider the excitations on an anisotropic ground state of a 
condensate in a cylindrically symmetric harmonic trap with trapping potential 

W) = l -muo 2 [x 2 + y 2 + (Xz)% (4) 

where A = uo z /uj r is the anisotropy parameter of the trap. Our numerical solutions 
for the ground state of the GPE and for the corresponding BdG equations (|3|) are 
expressed in units of the harmonic oscillator length r = \J%j 2mu r and energy fko T . 
We will illustrate our discussions using the case of uj t = 2n ■ 75 Hz and non-linearity 
parameter C = NU Q /hu T rQ = 332 (unless otherwise stated) though we have tested our 
results for a wide range of C values of up to 300000. Because the wave function Vi 
is in general fairly much a mirror image of Ui, and \ui\ 3> (apart from the lowest 
excitations) ||, we present only the wave function Uj. 



3. Symmetries 

For isotropic harmonic traps the ground state of the GPE has I = 0, and so the GPE is 
completely separable and reduces to a one-dimensional problem in the radial coordinate 
r. The BdG equations (for excitations on the ground state) are therefore also spherically 
symmetric, and thus the total angular momentum and the z-component of the angular 
momentum are conserved. This allows a separation of the angular dependence in the 
usual fashion by writing the quasiparticle amplitudes as u(r) = u T {r)Yi m { , &,^p) and the 
equivalent for v(r), where YJm($, V 9 ) are spherical harmonics. Since the radial equation 
for the excitations does not contain any dependence on m, the solutions with the same 
I and different m are degenerate. 

For cylindrically symmetric harmonic traps a complete separation of variables is 
not possible. Although the ordinary Schrodinger equation is separable in cartesian or 



cylindrical coordinates [30], this separation is not possible for the GPE due to the non- 
linear term. Nevertheless, solutions of the form 

iP(r) = tp(p, z)e [m ^ (5) 

can be found, where p, ip and z denote the usual cylindrical coordinates. Then, neither 
the trapping potential nor |-?/>(r)| 2 in the non-linear term of the GPE depend on the 
azimuthal angle ip and thus L z and the operator Ho + NUo\ip(r)\ 2 of the GPE commute, 
so that m c is a good quantum number. Since no further separation is possible, the 
equation must be solved in the two variables p and z. The ground state solution of the 
GPE equation corresponds to m c = 0. 

Correspondingly the normal modes of the BdG equations will also have specific 
angular momentum compositions if ip(r) is given by (|5|). If Ui(r) is an eigenfunction of 
L z with eigenvalue hm, then i>j(r) will be an eigenfunction with eigenvalue h(m — 2m c ) 
|3~T|. Excitations on the ground state (m c = 0) with ±|m| are degenerate because m 



enters quadratically into the BdG equations in cylindrical coordinates. 

The axially symmetric trap potential also has a reflection symmetry with respect 
to the x-y plane, and thus the solutions to the BdG equations can be chosen to have a 
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well-defined parity ||. In the isotropic case, where / and m are good quantum numbers 
for the excitations, the parity is simply given by II = (— l) l ~ m . Hutchinson and Zaremba 
showed that, as the trap geometry is altered to an anisotropic one, the degenerate modes 
split into branches with different excitation frequencies, which emerge continously from 
the isotropic case (see figure 1 in reference This allowed them to associate a number 
I with each branch depending on where it originates in the isotropic limit, even though 
total angular momentum is no longer a good quantum number. 



4. Classification of Excitations on the Ground State 



4-.1. Thomas-Fermi limit 

In the Thomas-Fermi limit a complete separation of variables is possible. Stringari 



TT| has shown that by casting the GPE into hydrodynamic form and then taking 
the Thomas-Fermi limit in the hydrodynamic regime (i.e. energies much smaller than 
the chemical potential) a second order wave equation can be derived to describe the 
elementary excitations. Fliesser et al |T1J recognized some underlying symmetries 
of this equation by identifying three operators which commute with each other, and 
introduced three corresponding quantum numbers (n, j, m) that classify the solutions 
completely. An explicit separation of the wave equation was achieved in cylindrical 
elliptical coordinates £, n and <p, and in terms of these variables the quantum numbers 
represent 

n: order of polynomial in £ and n 

j: index to label different eigenvalues for fixed n and \m\ ; 

j runs from to iV = 1 + int | 
m: ^-component of angular momentum. 

The Thomas-Fermi approximation is known to be very accurate as a limiting case of 
the full solution if Na/r ^> 1, where r is the harmonic oscillator length. In particular, 
the shape of the ground state wave function is very well approximated. Although the 
solutions of the full BdG equations do not strictly conserve these quantum numbers, we 
find that they exhibit in general the same patterns and symmetries, and we will show, 
in the appropriate regime, how the family classification scheme we develop in this paper 
can be related to n,j, m. 



4-2. Families in the isotropic case 

We first consider the isotropic case because then the patterns of the different mode 
families can be easily described in terms of Legendre polynomials. As the trap geometry 
is changed from spherical to cylindrical symmetry these patterns are continously 
modified, being squeezed in the direction of the stronger confinement, but the basic 
character remains recognizable. 
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In cylindrical coordinates the solutions for the excitations in an isotropic trap, 
u(r) = u r (r)Yi m ($, <p), can also be separated as u{r) = u(p, z)e im</3 . The relation 



Y lm {$,ip) 



2l + Hl - m) \ Plm (cos^ (6) 



/2^\ 2 (l + m)\ 

between the spherical harmonics Yimf^, ip) and the Legendre polynomials P^cos'i?) 
shows that u(p, z) is essentially the radial function u r {r) modulated by the Legendre 
polynomial P; m (cosi9), where cosi? = z(p 2 + z 2 )~ l l 2 in cylindrical coordinates. In fact, 
it turns out that the general shape of the families is determined by the symmetries of 
the Legendre polynomials. We now show that the family classification suggested by 
Hutchinson and Zaremba can be generalized, in the isotropic CclSC, clS follows. First we 
assign a principal family number which is given by 

f = l-\m\ + l, (7) 

and then an additional number characterizing the radial function is needed to complete 
the classification into families. We shall introduce the nodal family number n r for this 
purpose, which in the isotropic case is simply the number of nodes in the radial function. 
The family is given by the pair (f,n r ), which together with the magnetic quantum 



number m uniquely specifies any mode. In section [4.3| we will show how this family 
classification generalizes to the anisotropic case. 

We illustrate the spatial character of the family assignment by considering first the 
excitation modes with no radial node (n r = 0). We begin with the case m ^ 0, which 
we illustrate in figure [3] with contour plots of full numerical solutions of u(p, z) for the 
specific case of the degenerate modes of the lowest I = 3 excitation with m = 3,2, 1. 
Their principal family numbers are / = 1, 2, 3 respectively. The cylindrical symmetry 
means that the p-z dependence can be found on any plane through the z-axis (we have 




Figure 1. General shape of mode families 1 to 3 with no radial node. Contour plots in 
the x-z plane of the quasiparticle amplitude u(p, z) are given for the degenerate modes 
I = 3, (a) m = 3, (b) m = 2, (c) m = 1 modes. 
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chosen the x-z plane), and the modes are symmetric with respect to the x-coordinate. 
Since the radial function is the same for each of these modes, the relative overall shape 
is determined by the Legendre polynomials. The important property of the Legendre 
polynomials P; m (cosi?) for our purposes is that they have n# = I — \m\ nodes between 
< d < n. Thus, the number of angular nodal surfaces n$ beween < < ir, i.e. 
surfaces of zero density that are characterized by a constant value of $ in the isotropic 
case, determines the principal family number / since / is given by equation (0) as 
/ = I — \m\ + 1 = n# + 1. We note also that for m ^ all Legendre polynomials are zero 
along the z-axis, but this is not a nodal surface. Because the sign of the wave function 
changes as it crosses a nodal surface, family 1 members have even parity, family 2 have 
odd parity, and in general the parity of the mode is related to the principal family 
number by II = (— l)^ -1 . 

The m = member of each family has a shape that derives from the Pio Legendre 
polynomial. We call it the anomalous member of the family, since its shape differs from 
other members of the family only in that it is non-zero along the symmetry axis, which 
does not change the character of the excitation significantly. We illustrate the shape of 
the anomalous modes of the families (2, 0) and (3, 0) in figure |2|. We note also that the 

Family (2,0) Family (3,0) 




x-position x-position 



Figure 2. Contour plots in the x-z plane of the amplitude u(p, z) of the anomalous 
first members of family 2 and 3 (m = 0). 

anomalous member of family (1, 0) is the ground state, which is a solution of the BdG 



equations p| . 

The case where the radial function has a non-zero number of nodes (i.e. n r ^ 0) 
can now be easily visualized. The principal family number / determines the number 
of angular nodal surfaces (/ — 1) between < d < tt, while n r determines the number 
of radial nodal surfaces, which intersect the angular nodal surfaces. In the isotropic 
case they are spherical and centered on the origin. In figure |3| we illustrate the first two 
modes having one node in the radial function, which both belong to the family (1, 1). 
The mode in figure ^| (a) is the anomalous first member of this family (m — 0), while 
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Figure 3. Family 1 with one radial node (n r = 1). (a) Anomalous first member (I = 0, 
m = 0), (b) general shape (here I — 1, m = 1). 

all other modes of this family have the general shape shown in figure |3| (b), which can 
be recognized as the same shape as in figure p] (a), but with one radial node. 

We stress that all members of the same family, apart from the anomalous one, 
have the same general shape, i.e. the same number of peaks (in the contour plot) in 
similar spatial distribution. The main qualitative difference between modes of the same 
family is that the peaks move radially outwards and become narrower in both radial and 
azimuthal direction with increasing eigenfrequency. We also note that in the isotropic 
case the principal family number / can be readily obtained from the contour plots by 
counting the peaks in the half-plane x > if there is no radial node (n r = 0), or counting 
the peaks in the region that is enclosed by the z-axis and the n T — 1 radial nodal surface 
if n r > 0. For m = all peaks centred on the z- axis must also be included. 

To illustrate the mode classification by family and m value, we list in table [I] the 
18 lowest lying modes for the C = 332 case in an isotropic harmonic trap. 

4-3. Anisotropic case 

The main value of the concept of families is in its extension to the anisotropic 
cylindrically symmetric case. Hutchinson and Zaremba identified the first four families 
by the dependence of the eigenvalue on trap anisotropy. Here, we show that the mode 
topology determines the family. 

In figure |] we have plotted the quasiparticle wave functions for three families: the 
anomalous member of the (3, 0) family, and the (2, 1) and (3, 1) families, for the case 
of prolate, spherical and oblate traps. These graphs illustrate the general features that 
are found in the anisotropic case of cylindrical symmetry. We see that just as for the 
ground state of the GPE the quasiparticles are squeezed in the direction of the stronger 
confinement of the trap and expand in the other direction. For A < 1 (prolate) they 
are compressed in the p-direction and for A > 1 (oblate) they are squeezed in the z- 
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Table 1. Lowest quasiparticle modes of a condensate in an isotropic trap for C = 332, 
listed by family. 



Mode Family Mode Family 
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direction. This distortion has no effect on the zero line along the symmetry axis, but the 
nodal surfaces of the isotropic case are distorted. The radial nodal surfaces are no longer 
spherical, but are changed to a shape that approximately follows the equipotentials of 
the trap. In all cases the total number of nodal crossings along the positive half of the 
axis of strong confinement remains exactly n r . The angular nodal surfaces can no longer 
be parameterized by constant and they do not meet at the origin: in the prolate case 
they generally intersect the symmetry axis at various (non-zero) values of z, while in the 
oblate case they generally intersect the z = plane in circles of different radii. We will 
call the former planar nodal surfaces and the latter cylindrical nodal surfaces, which is 
a good description in highly anisotropic traps as we will see in section |5.1| . Exceptions 
are rare and occur for the case where the trap is close to spherical. Then, it is possible in 
the prolate case that adjacent planar nodal surfaces join just before meeting the z-axis, 
while in the oblate case it is possible that adjacent cylindrical nodal surfaces join just 
before the z = plane. 

The nodal crossings of the isotropic case, which are at the centre of four peaks 
of alternating sign, can also change character in the anisotropic CclSG, clS illustrated in 
the right hand column of figure |] (a). There we see that pairs of peaks of the same 
sign can begin fusing; in other words the crossing has become a saddle point, or ' anti- 
crossing 1 . Despite these distortions, the character of the isotropic excitations remains 
clearly evident in the anisotropic case, and thus a family assignment can be made. The 
first step is to determine the number of radial nodes (by inspection along the positive 
part of the axis of strong confinement), which gives the radial family number n r . Next, 
if there is no radial node (i.e. n r = 0), we count the number n p of distinct peaks in a 
single quadrant (including peaks that are centered on the axes). If a radial node does 
exist, we identify the first radial nodal surface, and count the number n p of distinct 
peaks in a single quadrant inside this radial nodal surface. In either case, the principal 
family number is then given by / = 2n p — 1 for even modes and / = 2n p for odd modes. 
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Family (3,0) Family (2,1) Family (3,1) 




-8 1 \ \ 

-8 -4 4 8 -8 -4 4 8 -8 -4 4 8 
x-position x-position x-position 

Figure 4. Effects of trap anisotropy on quasiparticle shapes. Left hand column / = 3, 
n r = 0, to = mode; middle column / = 2, n r = 1, m = 1 mode; right hand column 
/ = 3, n T = 1, m = 1 mode. The anisotropy parameter is (a) A = 0.5, (b) A = 1, (c) 
A = y/8. 

The n r = 1 nodal surface can be traced out by following a path from the initial point 
on the axis of strong confinement, through successive crossings or anti-crossings, until 
the axis of weak confinement is reached. When anti-crossings are encountered, the path 
continues over the saddle (along the line of minimum amplitude) to the nodal line leaving 
the anti-crossing opposite to the entry. Of course this line (or radial surface) is no longer 
strictly nodal, but nevertheless serves to determine a region which is characteristic of 
the family. Just as for the isotropic mode is uniquely identified when the family 
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assignment (f,n T ) is given together with the magnetic quantum number m, and the 
parity is given by IT = (— 

5. Ordering of Quasiparticle Eigenfrequencies 

In this section we show how the family classification can be related to the energy 
ordering of the quasiparticles for given m, and we provide a simple model to explain the 
relationship. 

5.1. Full solutions of Bogoliubov-de Gennes equations 

In figure |5] we illustrate the shape of the first few family members with n r = and 1, 
for the case of very prolate and very oblate traps (figures [| (a) and (b) respectively). 
We have chosen to represent the m = modes, because the mode shapes are slightly 




-6 6 -6 6 -6 6 -6 6 _4 q 4 

x X 



Figure 5. Highly anisotropic quasiparticle eigenstates in (a) prolate trap with 1/A = 7 
and (b) oblate trap with A = 7. The cases shown are the / = 1,2,3,4 anomalous 
(m = 0) family members with n r = and 1, respectively. The contour plots show 
contour lines only at ±0.001. 

simpler than the m^fl modes (which differ only in having a zero along the z-axis). In 
highly anisotropic traps, the family shapes follow very well defined patterns, as we can 
see, and it is convenient to discuss the prolate and oblate cases separately. 

Prolate case: For the prolate case an increase in the principal family number / simply 
adds an additional planar nodal surface perpendicular to the z-axis. Changing n r from 
to 1 adds a radial nodal surface which appears in this projection as two lines symmetric 
about, and almost parallel to the z-axis. The energy ordering of these modes is plotted in 



Elementary excitation families and their frequency ordering in cylindrical BECs 1 1 

figure ^ (a) where we see that, for a given anisotropy, the energy ordering initially follows 
the family assignment (/, 0) with alternating even and odd modes. This sequence is 




Figure 6. Frequency spectrum versus anisotropy of the low- lying n r = modes and 
the first two n r = 1 modes (m = 0) in (a) prolate, (b) oblate geometry. The even and 
odd modes with n r — are shown as dotted and dashed lines, respectively, and the 
n r = 1 modes as dashed-dotted lines. The solid line shows the condensate chemical 
potential \i. The mode eigenfrequencies are measured relative to the condensate 
eigenvalue fi. 

interrupted by the (1, 1) mode, at an energy near 2, where it has become more favourable 
to have a single radial nodal surface than many planar ones across the narrow dimension. 
We note that the higher the anisotropy (i.e. the larger 1/A) the more (/, 0) modes fit in 
before the (1, 1) mode becomes favourable. 

Oblate case: Figure |5| (b) shows that in the oblate case it is natural to separate the 
sequences of even and odd modes. Then, successive members of each sequence (for fixed 
n r ) are obtained by adding an additional cylindrical nodal surface about the z-axis, 
which in this projection appears as a pair of nodal lines parallel and symmetrically 
displaced with respect to the £-axis. Changing n r from to 1 adds again a radial nodal 
surface, which appears here as a pair of nodal lines parallel to and symmetrically placed 
about the z = plane. 

In figure |6| (b), where the energies of the oblate modes are plotted as a function of 
the anisotropy A, the distinction between the even and odd modes is clearly revealed. 
The odd modes are shifted as a group to higher energies than the even modes (reflecting 
the energy cost of the nodal surface through z = 0), so that the energy ordering at a 
given anisotropy no longer simply alternates between even and odd. We see too that 
the even modes have essentially constant spacing except for very low A, and the spacing 
of the odd modes, while compressed for low excitation numbers, becomes equal to the 



Elementary excitation families and their frequency ordering in cylindrical BECs 12 



spacing of the even modes for higher excitation numbers. Once again, as for the prolate 
case, the higher the anisotropy the more (/, 0) modes fit in before the first n T — 1 mode. 

5.2. Comparison with harmonic oscillator solutions 

For higher excitations the quasiparticle amplitude Vi is negligible compared to Ui and 
the BdG equations (H) reduce to the eigenvalue problem 



which has the form of a single particle equation 0, |12| . The presence of the condensate 
gives an effective potential that is broadly harmonic, but with a repulsive dimple in the 
middle, of height approximately ji. At high excitations, the energy levels for equation (|8|) 
will be essentially those of the harmonic oscillator, while for mode energies comparable 
or less than /i, the presence of the dimple will shift the energy levels upwards, effectively 
compressing them. 

If the condensate density is ignored in equation (||) the equation is simply that of 
the harmonic oscillator and is separable in cylindrical coordinates. The energy levels 
are given in units of huj r by 



where n p and n z are the number of nodes in the p- and z-wave function respectively 
{rip = 0,1,..., n z = 0,1,...). Equations @ and (|T0|) provide considerable insight into 
the behaviour of the solutions we have obtained from the full BdG equations. 

In the prolate case, increasing n z corresponds to increasing / in the BdG solutions, 
while increasing n p corresponds to increasing n r in the BdG solutions. For A < 1 it 
becomes energetically favourable to increase n z rather than n p , and we can fit roughly 
1/A nodes in the ^-direction before it is energetically favourable to put a node in the 
p-direction. In figure [7] we show from calculations of the BdG equations the number 
jyufi) f (y^ modes of lower energy than the first n r = 1 mode and compare this 
qualitatively to the predictions from equations (J|) and (llOf ). We see that while there 
is relatively good agreement at C = 332, the full BdG calculation predicts more (/, 0) 
modes to fit in, and the discrepancy increases at higher values of C. This is due to the 
fact that in the full equations the presence of the condensate causes the spacing of the 
lower levels (up to an energy of order \i above the condensate) to be compressed, and \x 
increases with C. 

In the oblate case, the agreement between the harmonic oscillator solutions and 
the BdG solutions is much better. As explained in the previous section, it is useful to 
separate the even modes (where n z corresponds to 2n r , and n p to (/ — l)/2) and the 
odd modes (where n z corresponds to 2n r + 1, and n p to (/ — 2)/2). Since A > 1, it 
is now clear from equations (^) and ( |T0| ) that increasing n z , i.e. increasing n r in the 
BdG solutions, is more energetically costly than increasing n p , which corresponds to 
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(9) 
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Figure 7. Comparison of energy ordering predictions of full BdG and harmonic 
oscillator solutions in a prolate trap. Plotted points show N^'°\ the number of n I = 
modes with energy lower than the (1, 1) mode. All modes are m = 0. 



increasing /. Figure || provides a quantitative comparison of the predictions of the 
harmonic oscillator solutions and the BdG solutions for the number of (/, 0) modes of 
lower energy than the first n r = 1 mode. We have separated the even and odd modes 
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Figure 8. Comparison of energy ordering predictions of full BdG and harmonic 
oscillator solutions in the oblate case for (a) the number A^iven of even n r = modes 
with energy lower than (1, 1) mode, (b) the number of odd modes with energy 

lower than the (2,1) mode. All modes are m = 0. 

so that figure § (a) shows the number of even modes with eigenfrequency lower 

than that of the (1,1) mode, while figure [8] (b) shows the number Nf£jp of odd modes 
with eigenfrequency lower than that of the (2,1) mode. The agreement is excellent for 
the C = 332 case and still very good for the C = 10000 case and is much better than 
for the prolate case of corresponding asymmetry, because the chemical potential \x lies 
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lower relative to the eigenfrequencies in question, as can be seen in figure |6|. 

Finally, we remark that the magnetic quantum number m simply gives an energy 
offset in equations (||) and fllDD, i.e. modes with a given number of nodes have higher 
energies for higher m. It is clear therefore, that the results discussed above for m = 
will also apply for m/0. The energy offset by m also explains that for a given family 
(/, n r ) (i.e. a given number of nodes) the energy increases with increasing m (see table 



6. Discussion 



We have shown that the concept of mode families introduced by Hutchinson and 
Zaremba || can be systematically defined and extended to include all excitation modes 
of cylindrically symmetric anisotropic traps. The family assignment (/, n r ) determines 
the topology of any mode with m/0, and the m = mode differs only in being non- 
zero on the symmetry axis. In the regime where the treatment of Fliesser et al is 
valid, we can relate their quantum numbers n, j to our family assignment numbers / 
and n r as follows: 

f = n-2j + l, n T =j. (11) 

The similarity of shape of modes of the same family explains the similarity of 
frequency dependence on anisotropy for modes in a given family found by Hutchinson 
and Zaremba; changes in trap geometry affect every mode in the family in much the 
same way. We have also shown how the energy ordering of the quasiparticle excitations 
is related to their family shape and provided a simple model that explains the behaviour 
in terms of harmonic oscillator eigenstates. 
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